Resolved atomic interaction sidebands in an optical clock transition 
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We report the observation of resolved atomic interaction sidebands (ISB) in the ^'^Sr optical clock transition 
when atoms at microkelvin temperatures are confined in a two-dimensional (2D) optical lattice. The ISB are a 
manifestation of the strong interactions that occur between atoms confined in a quasi-one-dimensional geometry 
and disappear when the confinement is relaxed along one dimension. The emergence of ISB is linked to the 
recently observed suppression of coUisional frequency shifts in |T|. At the current temperatures, the ISB can 
be resolved but are broad. At lower temperatures, ISB are predicted to be substantially narrower and usable as 
powerful spectroscopic tools in strongly interacting alkaline-earth gases. 



Experimental efforts in control of ultracold alkaline-earth 
atoms (AEA), such as Sr or Yb, have led to remarkable devel- 
opments in optical atomic clocks yj|2J, which are approach- 
ing the accuracy of single ion standards O. Fermionic AEA 
are also beginning to attract considerable attention in the con- 
text of quantum information processing |4| and quantum sim- 
ulation |l5]|6l. Many of these applications require reaching 
the strongly interacting regime, in which interactions domi- 
nate over all relevant dynamical energy scales. In this Letter, 
we demonstrate the onset of strong interactions in a low di- 
mensional gas of '^^Sr atoms by observing resolved interaction 
sidebands (ISB) in Rabi spectroscopy of the ^5*0 — '^Pq optical 
clock transition. 

Rabi spectroscopy has served as an excellent probe for in- 
teracting alkali gases. For example, it was used to detect and 
characterize Rose Einstein condensation in spin polarized hy- 
drogen |7|, to directly measure the fermion pair wave func- 
tion in the BEC-BCS crossover fSl, and to resolve the Mott 
insulator shell structure in a bosonic lattice system |9|. Here, 
we show that the ISB structure that emerges during optical 
Rabi interrogation can become a powerful spectroscopy tool 
to uniquely probe lattice-trapped AEA. At microkelvin tem- 
peratures and loading conditions of mainly one or two atoms 
per lattice site, we observe a single ISB, which can be resolved 
from the carrier but is broad. However, in the quantum degen- 
erate regime, the ISB are expected to become substantially 
narrower and could thus lead to a precise determination of 
the s-wave ^Sq — ^Pq interaction parameters (currently, the 
— ^Pq scattering lengths for fermionic AEA remain un- 
known). In addition, we show that the variation in ISB spectra 
between sites with different numbers of atoms allows for a 
precise characterization of atom number occupations in inter- 
acting AEA systems. 

When the ^Sq and ^Pq electronic degrees of freedom 
within a single fermionic atom are represented as an effective 
pseudo-spin 1/2 (g and e), two nuclear-spin polarized atoms 
form a simple model for understanding the appearance of ISB 
(see Fig. [T]). We consider atoms trapped in a tube geome- 
try with a weak trapping frequency uz along Z and a strong 
transverse confinement hojx.Y ^ kBTx,Y, large enough to 
forbid any transverse dynamics [Tx{Ty) is the temperature 
along X{Y)]. We assume the atoms are initially prepared in 
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FIG. 1: (Color online) Schematic diagram of Rabi spectroscopy 
in the strongly interacting regime. When 5 ~ 0, only the non- 
interacting triplet states (S = 1) are accessible, giving rise to an 
interaction-free carrier However, when S ~ Un, the singlet |0, 0) be- 
comes resonant, and a slightly inhomogeneous coupling (AQ^ ^ 0) 
can induce population transfer to it, as manifested in an ISB. 



the same pseudo-spin state {\gg)) and interrogated by a lin- 
early polarized laser beam with bare Rabi frequency fl^ and 
detuning from the atomic resonance S. The Pauli exclusion 
principle forces atoms in identical internal states to occupy 
different vibrational levels along Z, n = (ni, 712), since their 
spatial wave function must be antisymmetric. Being identical 
fermions, the atoms initially experience zero s-wave interac- 
tion. If the atoms are coherently driven to e (i.e. if they expe- 
rience the same Rabi frequency Clu = (f2„j + f2„2)/2), then 
the effective spin must remain symmetric under exchange. 
Consequently, during the excitation the atoms will not ex- 
perience any s-wave interactions. We refer to states with 
spin-symmetric wave functions, \gg) = \S ~ 1,M ~ — 1), 
|ee) = |1,1) and {\eg) + \ge))/V2 = |1,0), as the triplet 
states, S = 1 ifTOlfTTI . 

However, in optical transitions, a small component of the 
probe beam along Z leads to a slightly different Rabi fre- 
quency VLn [rfz) ^^'-h mode, with rjz = fczaf^/v^ the 
Lamb-Dicke parameter, m the atom mass, kz the component 

of the probe laser wave vector along Z and a^^ = \J^;^^ the 

Z harmonic oscillator length. If Ari,i = (il„j — ri„2)/\/2 
is not zero, then the optical excitation-induced inhomogeneity 
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can transfer some of the atoms to the antisymmetric pseudo- 
spm smg let state, {\eg) - \ge))/V2= |0,0). The mteractmg 
singlet is separated from |1, 0) by an interaction energy Ufi. 

In the limit where 5 ~ fi^j <C J7s, a small A57f; cannot 
overcome the energy cost required to drive the transition to the 
singlet (see Fig. [T}. Evolution into the singlet is blocked and 
the atoms remain non-interacting. However, if we increase the 
detuning to 5 ^ Un, the transition |1, —1) — > |0, 0) becomes 
resonant. Here, a small Ail,-j can efficiently transfer the atoms 
to the singlet, giving rise to an ISB in the lineshape. Note that 
if instead of preparing all the atoms in g we prepare the atoms 
in e, the ISB occurs at 5 ^ — ?7,-j. 

The ISB resemble motional sidebands that develop in spec- 
troscopy of trapped ions when the trapping frequency is much 
larger than the recoil energy, the so called Lamb-Dicke regime 
IIT2I . In this case, the carrier becomes free of Doppler shifts 
up to the corresponding "Stark shift" 113], and all motional 
effects are manifested in the sidebands. In our case, the ISB 
are pushed away from the carrier due to the large energy sep- 
aration between interacting and non-interacting states, leav- 
ing the carrier free from interaction effects, up to a correc- 
tion proportional to Arti/Un. Recently we have observed the 
suppression of coUisional frequency shifts to the 10~^^ level 
under similar tight trapping conditions |[T] . The observation of 
ISB unambiguously confirms this suppression to be the con- 
sequence of reaching the strongly interacting regime. 

Our experiment employs nuclear spin-polarized (mj = 
+9/2) *^Sr atoms loaded in two lattice configurations: a one- 
dimensional (ID) vertical lattice which creates an array of 
quasi-two-dimensional traps (pancakes) along Y with weak 
confinement along the X and Z directions, and a deep 2D lat- 
tice which creates an array of quasi- ID traps (tubes) in the 
X-Y plane and has weak confinement along Z. To prepare 
the atomic system, we laser cool ®^Sr atoms to about 2 /iK 
inside a magneto-optic trap based on the weak ^5*0 — ^Pi 
transition and then load them into a ID vertical lattice. The 
spatial distribution of occupied ID lattice sites is determined 
by the vertical extent of the MOT cloud, which is approxi- 
mately Gaussian with a standard deviation ay — 30 /im. For 
the 2D lattice confinement we then adiabatically ramp up the 
horizontal lattice (along X). To remove any atoms trapped in 
the ID vertical lattice outside the 2D intersection region, we 
ramp the vertical lattice off and then back on. The number 
of horizontal lattice sites occupied is determined by the radial 
temperature of the vertical lattice. Atoms load into ^100 rows 
of tubes uniformly distributed along Y, while the columns dis- 
tributed along X are loaded according to a Gaussian distribu- 
tion with standard deviation an of 6- 10 /im. After forming the 
lattices, we perform Doppler and sideband cooling using the 
^Sq — ^Pi, F — 11/2 transition. Simultaneously, atoms are 
optically pumped to the m/ ~ +9/2 ground state magnetic 
sublevel, using CT+-polarized light on the — "^Pi, F = 9/2 
transition, propagating along a bias magnetic field parallel to 
Z. We perform spectroscopy of the ^Sq — ^Po transition using 
a narrow linewidth laser propagating along Y. 

We quantify the number of atoms by detecting fluorescence 



on the strong ^5*0 — '^Pi transition at 461 nm. Temperature 
information and trap frequencies are extracted from Doppler 
spectroscopy along Z and vibrational sideband spectra fTlfT?). 
We estimate that we loaded ^7000 atoms in the ID lattice 
distributed in ^420 pancakes and ^17 atoms per pancake. We 
determine uy = 27r x 80 kHz and ujx,z = Itt x 500 Hz, 
and Tx^Y,z ^ 4 ^K. For the 2D lattice, we loaded --3000 
atoms and estimate that 30% of the populated tubes are doubly 
occupied. We determine ujx = 27r x 110 kHz, = 27r x 70 
kHz and ujz — 2tt x 800 Hz at trap center, and Tx.y.z ~ 4.5 

ISB develop in the lineshape only in the parameter regime 
where 7 = Eint/{^)T » 1. Here, Ei,,t = {N - 1)\{U)t\/2 
is the mean interaction energy per particle and N repre- 
sents the number of atoms per lattice site. {U)t character- 
ize the strength of two-body interactions and depends both on 
the confinement volume and temperature as |17|: {U)t ~ 

with u = 4a-^^- 
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the singlet ^6*0 — '^Pq scattering length (|6l and 1) and 
-d functions characterizing the temperature dependence, i.e. 
-dix > 1) - 2^{x > 1) 1 and < 1) - d{x < 
1) — ?► ^/x. We used (0)t to denote the thermal average of 
O. In our experiment, the — ^Pq spectroscopy is gener- 
ally performed by using a 80-ms pulse and the "7r-pulse" laser 
intensity is set to achieve maximum excitation fraction. We 
estimate (fi)T ^ 27r x 6.25 Hz. 




FIG. 2: (Color online) Measured Rabi lineshapes for 17 to e inter- 
rogation (normalized with respect to peak height). The main Panel 
shows lineshapes taken in a ID lattice (black triangles) and in a 2D 
lattice (connected green circles) when the transition is excited with 
a TT-pulse with 80 ms duration. The noninteracting, homogeneous 
Rabi lineshape for an 80 ms 7r-pulse is also shown (orange line). The 
inset shows lineshapes measured in a 2D lattice using a 160 ms 2n- 
pulse (connected red diamonds) and a 160 ms 7r-pulse (connected 
blue squares). Each lineshape was measured under the trapping and 
temperature conditions described in the main text. The data from 
each lineshape is collected into bins 1-5 Hz wide. 

Using the experimental trapping conditions and tempera- 
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tures and assuming a moderate a^g = — 70ao (uq the Bohr 
radius), we obtain ^ 10 in the doubly occupied sites at 
the center of the 2D lattice. On the contrary, for the central 
pancake of the ID lattice with a mean number of 17 atoms 
and the same scattering length, 7 is reduced to 7^^ ^ 0.6. 
Based on these values, one predicts the development of ISB 
only in the 2D lattice geometry, a prediction confirmed by our 
measurements. We measured a series of lineshapes in the ID 
and 2D lattices. In the 2D lattice case we also varied the pulse 
area and Rabi frequency. The resulting data are shown in Fig. 
|2] To remove some of the statistical noise we performed mul- 
tiple scans across the transition resonance under each experi- 
mental condition. For each scan, the drift of our ultra stable 
laser was canceled to <50 mHz per second, and the laser was 
stepped by 2 Hz each experimental cycle (^1.5 sec). The cen- 
ter of each scan was determined by a fit to a Lorentzian. Each 
data set consists of ^20 concatenated scans and the direction 
of the scan was alternated to reduce systematic effects due to 
residual drift. We separated the data into bins and calculated 
the mean and standard error of the mean for each bin. At our 
4 — 5 /iK temperatures, the ISB in the 2D lattice are notice- 
able but nevertheless broad. To improve resolution, we reflect 
the binned lineshape at positive detuning (free of ISB when 
fljTg < 0) and subtract it from the corresponding bin at nega- 
tive detuning. The resultant lineshapes with a cleaner ISB are 
shown in Fig. [4] 

The presence of ISB in the tightly confined geometry can 
be understood quantitatively by considering an isolated tube 
with N atoms. We characterize the tube dynamics by defin- 
ing a set of effective spin operators, S'^'^'^, in the {e, g} ba- 
sis |TT|. Here the subscript rij is drawn from a fixed set 
ft = {rii, 712 . . . riAf} of initially populated vibrational modes 
along Z. The description of the system in terms of effec- 
tive spin operators is valid provided those initially populated 
modes remain singly occupied by either a g or an e atom dur- 
ing the excitation process. In the rotating frame, the Hamilto- 
nian of the system up to constant terms becomes IITTI . 
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— X^jLi ^nf collective spin operators. The quantity 



Un,,n^, = ud{^^)d{-^)In,,n^, measures the Strength 
of the interactions t5etween two atoms in the antisymmetric 
electronic state. /„^ ,„^, is an overlap integral between har- 
monic oscillator modes along Z ifTTI . Only the Z— mode dis- 
tribution is treated exactly since at current temperatures only 
a few transverse excited modes are populated. The population 
of transverse modes is accounted for as a renormalization of 
the interaction parameter 

The interaction part of the Hamiltonian is diagonal in the 
collective spin basis |5, 7\f), 5* = 0(i), . . . , N/2 and |A/| < 
S. For N — 2, the spin basis is spanned by the triplet states 
|1, M) and the singlet |0, 0). Among the collective states only 
the S = N/2 states are noninteracting. States with S < N/2 
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FIG. 3: (Color online) Excitation fraction vs. 5 at T = for lattice 
sites with N = 2 - 5 atoms. The ISB reveal the 5" = N/2 - 1 
spectrum and can be used to probe occupation number in a lattice. 
The parameters used for the plot were rjz = 0.4, fl^ = 27r x 5 Hz, 



t = 1.57r/n" and u = 2tt x 2800 Hz. 



experience a finite interaction energy. In the presence of exci- 
tation inhomogeneity, S is no longer conserved and during 
excitation of the clock transition, atoms can be transferred 
mainly between S = N/2 and S = iV/2- 1 states. The N- 1 
coflective excitation modes \S = N/2 -1,M = N/2 -~l,q) 
(q ^ 1,...N - 1) with energies t/H can be accessed from 
the initially g-polarized state \N/2, —N/2) using S ^ Ui/^ . 
They give rise to — 1 sidebands in the lineshape (Fig.|3]l: 



with f{t,5,y) 
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2 J^j ^71, {N/2, N/2\S/j\N/2 - 1, N/2 - 1, g). The first term 
is associated with the carrier and the second with the ISB. For 
the case = 2, we have J/i'^ = Ufi and Afil'^ — Aflfi. 

At quantum degeneracy, only the lowest lying vibrational 
modes need to be considered n = {0, 1, ... — 1}. In this 
case, the interaction sidebands are very narrow. Assuming 
t = sir/n^, with s a numerical constant that describes the 
pulse duration, the peak width is just 2Vl^ / s. Measuring the 
interaction sidebands thus precisely determines a^g and char- 
acterizes the S = A^/2 — 1 spectrum, which has structure for 
N > 2. Since the peak height of the interaction sidebands 
scales as 77^, large probe inhomogeneities are required to ob- 
serve them at low temperature. 

In Ref. 1 9 1, interaction-dependent transition frequency 
shifts in RF spectroscopy were used to spectroscopically dis- 
tinguish sites with different occupation numbers, revealing 
the shell structure of the Mott insulator phase. The g and 
e states in AEA have no hyperfine structure, and hence RF 
spectroscopy is not applicable. However, as shown in Fig. [3] 
the optical analog of RF spectroscopy- the ISB spectroscopy 
- can distinguish sites with different atom numbers. This ca- 
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FIG. 4: (Color online) Finite temperature ISB in a 2D lattice (g to 
e): Symbols correspond to experimental data obtained after reflect- 
ing the binned lineshape at positive detuning (free of ISB) and sub- 
tracting it from the corresponding bin at negative detuning. Data con- 
ditions previously shown in Fig. [2] are displayed with the same sym- 
bols and colors. The data indicated by purple crosses corresponds to 
an 80 ms yr-pulse but with different trapping conditions (u}x,y,z = 
27r X {92, 60, 0.6} kHz at trap center, rx,y,z ~ {2.5, 2.5, 7.2} /xK). 
The inset shows the absence of ISB in a ID lattice. Solid curves cor- 
respond to theoretical predictions. The dip at S ~ —70 Hz in the 
measured ISB is due to residual mi — 7/2 atoms. 



pability can have important implications in quantum simula- 
tion [6 1 and information proposals [4| with AEA. 

At finite T, instead of a fixed set of populated modes, a 
thermal average should be taken. For simplicity and to cor- 
respond to our current experiment, we focus on the N ~ 2 
case. After evaluating the thermal average fTTI, we obtain 
the following approximate expression for the ISB, valid in the 
regime \6\ > fl^, ksTz > hujz and t = sn/fl^ < {^^)tI 
(These conditions are satisfied in our experiment): 
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Fig. |4] shows comparisons between the measured ISB and 
theoretical curves calculated using Eq. (|3]l. In the 2D lat- 
tice, sites near the wings of the lattice beams' Gaussian 
intensity profiles are significantly shallower than near the 
center of the beam intersection region, i.e lo±{X,Y) k, 
a;_L(0,0)Exp[-(X2 + r2)/yy2j^ ^j^j^ = ^ujxojy and 

W_L = y/WxWy ~ 30 nm the mean beam waist. This im- 
plies u varies from tube to tube; so to compare with exper- 
iment, it is important to average Eq. ([3| over the populated 
tube distribution. We assume a uniform distribution of dou- 
bly occupied tubes. With this procedure, we obtain lineshapes 
that reproduce the experimental data if we use a~g = — 280ao 
and r]z ^ 0.07(0.05) for the = 27r x 87(73) kHz cases. 
The value of a^Tg that fits the experimental data is outside the 
range of effective scattering lengths a'^* ^ —(35 — 50)ao re- 



ported in Ref. |1J. However the discrepancy is attributed to 
the fact that a~* in Ref UJ corresponded to the effective scat- 



tering length required to model the entire lattice array by a 
single tube at the trap center, and thus \a'^* \ should be smaller 
than |a^g|. In agreement with theory, the ISB's position and 
width are on the order of {U)t and their resonant peak heights 
(normalized to the carrier) decrease with decreasing fl^ and 
increase with pulse area stt. 



We determine the nuclear spin purity of the atomic sample 
to be > 97% by scanning the probe laser over the clock tran- 
sition frequencies for other nuclear spin states. The signature 
of the remaining 3% is the small resonant peak at (5 ^ 70 
Hz in Fig. [2] By comparing the amplitude of this small peak 
with the ISB, we conclude that the observed sidebands are not 
caused by s-wave collisions between unpolarized atoms. 

We also performed e to 5 Rabi interrogation by first trans- 
ferring all the atoms from 5 to e with a resonant pulse, then 
removing the remaining g atoms and finally interrogating the 
e atoms with an 80 ms tt -pulse. However, even in the tight 
geometry we did not see evidence of ISB in this case. The 
lack of ISB can be attributed to the absence of double occu- 
pancy in a single tube, possibly due to a combination of the 
following reasons. (1) A decrease in the number of interro- 
gated atoms due to the additional steps required to prepare the 
atoms in e, (2) A more favorable transfer from g to e in tubes 
with a single atom and (3) Inelastic two-body e — e losses due 
to ]3— wave collisions similar to the ones recently measured 
in KRb molecules ifTSi and in ®^Sr (s-wave) lfT6l . We mea- 
sured the e atom lifetime in a ID lattice to be 600ms. The 
measured loss rate, rescaled by the corresponding density in 
the tight confinement case, translates to decay rates about 30 
times faster, i.e. lifetimes ^ 20 ms. Given that there is a 10 
ms waiting time between the preparation and the clock inter- 
rogation, e — e losses seem to be the relevant mechanism that 
leads to the absence of ISB. 

.(3) 

The effect of tunneling between tubes could in principle 
modify the shape and height of the ISB (but not its mean posi- 
tion) and generate an additional frequency shift of the carrier 
The latter can be particularly relevant in the strongly interact- 
ing regime where collisional shifts would otherwise be sup- 
pressed. In order to investigate the effect of tunneling, we also 
took one set of data under conditions which reproduced those 
used in Ref. The presence of ISB in this lineshape, mea- 
sured with Tx.Y ~ 2.5 (see Fig.|4]i, completely rules out 
the possibility that these spectroscopic features are caused by 
tunneling since, under these conditions, we predict that tun- 
neling rates are suppressed by two orders of magnitude com- 
pared to Tx,Y ~ 4.5 iiK. 

The studies presented here demonstrate that tools devel- 
oped for clock experiments can provide precise understanding 
of strongly interacting many-body systems. This work was 
supported by ARO DARPA OLE, NIST, NSF, AFOSR, NRC, 
and Lee A. DuBridge and NDSEG fellowships. 
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Supplemental material 



Interaction parameters 



The (absolute value of the) mean interaction energy per par- 
ticle -Eint is {N—1)\{U)t\/% where is the number of atoms 
in the trap. 



We consider nuclear-spin polarized fermionic atoms with 
two different electronic degrees of freedom (e, g) confined in 
a harmonic trap with trapping frequencies ijJx,y,z and assume 
the atoms interact only via s-wave interactions. The matrix 
element that describes colUsions between two atoms {i — \,2) 
in the — {nix,niY, riiz) harmonic oscillator eigenmodes 
<I>1. {XWn„ (X)4>1, {Z) is given by 



(A.4) 



where R = {X,Y,Z), u = ^a'g^ 



h 



vV2"i+'»2 (ni!)(n2!) 



-, and Hn{z) are Hermite polynomi- 
als. 

At finite temperature, expectation values need to 
be calculated by averaging over all possible combi- 
nations of modes weighted according to their Boltz- 
mann factor: {U)t ^ u^i^m^M^) 



with i}{a) 



»("i+r.2) 



3(ni+n2) 



»(ni+"2) 



and i!}{a) 



V -c.(ni+n2) — ■ The i9 function is introduced 
since in our experiment atoms start polarized in g, and even 
at T = they are required to occupy different modes along Z 
to satisfy the Pauli exclusion principle. 

For J^j^ <C 1, eigenmodes with \niz — n2z\ ^ 1 domi- 
nate, so that 



K 



niz+n2Z 
\niz-n2z\ 



(A.5) 



where K is the complete elliptic integral of the first kind. Fur- 
thermore, we checked numerically that, for our parameters, 
ISB is not significantly perturbed by approximating K with 
K{0), so that Iniz,n2z ~ ^ /, ^ =r- Using this approx- 

\nr'^/\niz-n2z\ 

imation, one can show that '&{a <i; 1) ^ '!?(a <^1) ^ \foL. 



Interaction sidebands (ISB) at finite Temperature for two atoms 

(N=2) 

Here, we derive Eq. (3) in the main text by performing a 
thermal average over different axial modes. First, consider a 
given set of populated modes along Z,n = {n\,ni). If the 
atoms are in the strongly interacting regime, O3, the 

ISB is given by 



(A.6) 



with/(t,5,t/) = ^sin^(^^ 

At finite temperature, we need to evaluate {N'')t = 
J: \ e-^("i+"2) With a = ^ 

-^ni >n2 

{N^)t under the approximations Afi^j 



We evaluate 



nBT;g(rai-TO2) 
V2 



and 



To carry out the sum, we treat a{n\ — n-2) = x and 
a(ni + n-i) = y as continuous variables (valid if a ^ 

1). We also define k{x) = ^^^^^ = bx and U{x) = 

^. Then {N^)t « N^{x/a)e-''dx with N^{x/a) = 

sin^ [W^E±^EEM^^ . Noticing that ^ is 

sharply peaked at xq = we evaluate and k at 

x = Xq. We write x = Xq + x' and expand U in x' to first 
order. We then define p = \$,\x' / (26xq^^), so that 

91, 5/2 

(cVl+p2j, (A.7) 

where c = tbxo/2. For c <C 1, the integral in Eq. (S4) is 
approximately equal to ttc. This yields Eq. (3) in the main 
text. 



